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CURRENT DISTRIBUTION IN A THIN FILM SUPER- 


CONDUCTING STRIP TRANSMISSION LINE* 
by A. R. Sass^ and I. D. Skurnick 
Lewis Research Center 

SUMMARY 

In the long wavelength limit the current distribution in a thin film superconducting 
strip transmission line can be described by an inhomogeneous Fredholm integral equation 
of the second kind. When a fluxoid conservation derivation of this equation is considered, 
physical insight into the structure of the kernel follows naturally. An approximate ana- 
lytic solution to the integral equation is derived for a specified range of geometrical 
parameters commonly encountered in practice. The solution is obtained by making use of 
the Liouville-Neumann method of successive approximations, and by approximating the 
resulting series by a series involving powers of a defined coupling factor. It is shown 
that the critical current of the thin film superconducting strip transmission line, based on 
the calculations in this report and a critical current density hypothesis, is underestimated 
by less than 10 percent. 


INTRODUCTION 

Several authors have indicated that superconductive computer components, which are 
constructed in the form of thin film strip transmission lines, are advantageous from the 
standpoint of switching speed and miniaturization (refs. 1 to 6). Superconducting strip 
transmission lines are also useful in the transportation of electrical information within a 
superconducting computer due to inherent negligible loss characteristics and high group 
velocity (refs. 2 and 5). The latter is true only if the film thickness is larger than, or 
comparable to, the London penetration depth. It has also been shown that if the film 


*A condensed version of this report was published in J. Appl. Phys. , vol. 36, no. 7, 
July 1965, pp. 2260-2267. 

^Presently employed at the RCA Laboratories, Princeton, New Jersey. 
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thickness is less than the penetration depth the group velocity is appreciably decreased, 
which makes the strip line useful for delay line memory application (ref. 2). In all 
these devices it would be useful to be able to predict the total current which can be car- 
ried by the strip line before it becomes normally conducting. 

Several microscopic theories have been advanced which indicate that switching in a 
thin film is initiated by a critical current density (ref. 7). Cooper and Marcus have 
shown, independently, that the problem of thin film switching is complicated by the fact 
that the current density is not constant over the cross section of the film (refs. 8 to 10). 
By using a formal Green’s function approach and employing the London current-field re- 
lation, Cooper derived the general inhomogeneous Fredholm integral equation of the 
second kind describing the current density distribution in a single film. He also obtained 
a specialized equation for the case in which the film thickness is less than or equal to the 
penetration depth - the thin film case. Marcus derived an identical equation by combin- 
ing the Biot-Savart law and the London equation, and obtained computer solutions for the 
thin film case. These methods can be used to derive an integral equation for the cur- 
rent distribution in a strip transmission line. Due to the complexity of the kernel, 
analytic solutions to the integral equation have not been found to date for either the single 
film or strip line cases. 

It will be shown that when the integral equation is derived for the superconducting 
strip transmission line, by using the concept of fluxoid conservation and the London 
current-field relation, certain useful properties of the kernel can be readily deduced. 

The integral equation is solved by using the Liouville -Neumann method of successive ap- 
proximations. The properties of the kernel allow the approximate evaluation of the 
Liouville-Neumann series for a range of geometric parameters of practical interest. 

The analytic solution for the current density distribution in the strip line exhibits the 
same property of current peaking at the film edges that Cooper and Marcus found for the 
single film. The critical current of the system can be calculated from a knowledge of the 
current density in the strip line and the critical current density obtained from micro- 
scopic theory and/or independent experimentation. 


SYMBOLS 

The rationalized mks system of units is employed throughout the report, 
unit vector 
magnetic induction 
C constant determined by total current 
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thickness of film 
electric field intensity 
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magnetic field intensity 

magnetic field intensity at point i due to current element at point j 
total current 
current per unit area 
current per unit width 

kernel of one dimensional Fredholm integral equation 
£/d 

distance between centers of two films in a strip-line pair 

field point in x-y plane 

source point in x-y plane 

position vector denoting surface of conductor 

time 
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— y dimensionless field point in one-dimensional analysis 
w 

2 

- y’ dimensionless source point in one-dimensional analysis 
w 

width of film 


London penetration depth 


coupling factor = K(u, u')du’ 


(/3d) 2 (w/d)/8ff 

permeability of free space, 4 !tx 10" H/m 
electrical conductivity measured in mho/m 
magnetic susceptibility 


3 



STATEMENT OF PROBLEM 


It will be interesting to compare some general features of the current distribution in 
a strip transmission line for the cases where (1) the conductors are both normally con- 
ducting, (2) both conductors are superconducting, and (3) one conductor is supercon- 
ducting and the other is normally conducting. Consider the line structure in the form of 
two parallel cylindrical conductors shown in figure 1(a). For the sake of generality it is 
assumed, at first, that the cross section of each conductor is arbitrary. It is further 
assumed that (a) the wavelengths of the fields propagating along the structure are much 
larger than the transverse dimensions of the system so that a quasi-static analysis is 
valid, and (b) the current density is in the z-direction and therefore is only a function of 
x and y. The latter condition is a consequence of charge conservation for the quasi- 
static case (v • J ~ 0). 

It is convenient to first discuss the situation in which both cylinders are normally 
conducting. Faraday's induction law, expressed in integral form, is 




dS 


( 1 ) 



where ds - is an element of the contour 

bounding an area S. Since J = ctE in 

-5 

both conductors and x m « 1 ~ 10 
for most ordinary paramagnetic and 
diamagnetic materials, equation (1) 
can be rewritten as 

/ J • ds = - ap Q ± fj H • dS (2) 



(b) Enlarged view of conductor showing typical contour of integratioa 
Figure 1. - Array of infinitely long cylindrical superconductors. 


Applying Faraday’s law to the dotted 
contour in figure 1(b) and letting the 
path in the z-direction be of unit length 
yield 

J(r) - J(0) = c rp 0 ^ H(r").(£ z x d?’) 

(3) 

where r is the usual position vector 

in the x-y plane, £ is a unit vector in 
z 
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the z-direetion, and dr" is a differential vector line element of integration. The mag- 
netic field intensity at the point r" generated by current elements in both conductors is 
H(r"). For the static case, the right-hand side of equation (3) is zero and therefore the 
current distribution in each conductor is uniform. It is interesting to note that the dis- 
tribution in a normal conductor is uniform even if the neighboring conductor does not 
have a uniform distribution. An example of this point is the situation where the second 
conductor is a superconductor. This case will be described later in the report. 

For the case in which both conductors are perfectly conducting (cr = °°), the right- 
hand side of equation (3) is still zero in the static situation. However, at some time in 
the past a transient existed so that it appears that infinite current densities were gen- 
erated. This unrealistic situation is alleviated by stipulating that during the transient 
and afterwards the current flows in an infinitely thin layer at the perfect conductor sur- 
faces with a distribution such that H is time independent everywhere inside the body of 
the conductor. Starting from zero field initial conditions it is evident that H is zero 
everywhere inside of the cylinders. Thus the magnetic flux through any imaginary sur- 
face in the interior of either conductor is zero. It is easy to show that the magnetic field 
at a point F in a cross section, due to an infinitely long line current I, which intersects 
the plane of the cross section at r* is given by 



a z x (r - r») 



( 4 ) 


If there is to be no component of H(r) perpendicular to the surfaces, then at every point 
on the surfaces of both cylinders it must be true that 

5 ( ? sl,2) ' S l( r sl,2) =° 

where F g ^ and r g 2 are the position vectors denoting the surface of conductors 1 and 2, 
respectively, and a^r^ g) is the unit vector normal to the surface of the conductors. 

Since H^r s j 2 ) found by summing over the contribution to the field made by cur- 
rent elements at the surfaces of both cylinders, it follows that H vanishes everywhere 
inside of the cylinders if J (rjg), the current per unit width at the conductor surface, is 
the solution to the pair of coupled integral equations 
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dr sl 



= 0 



dr sl 




> 




(5) 


While these equations will not be solved here, two interesting consequences are apparent 
from their form alone. 

(1) The surface current distribution in conductor 1 is influenced by the cross- 
sectional shape of conductor 1, conductor 2, and the separation between them. 

(2) The surface current distribution in conductor 1 depends on the distribution in 
conductor 2. The distribution in conductor 1 is altered by the presence of conductor 2 
even if conductor 2 carries zero net current. 

The static case in which a = °° can be simulated by a dynamic situation in which the 
skin depth is much less than the transverse dimensions of the conductor. 

Consider the case in which both conductors are superconducting. In this situation it 
is no longer correct to use the expression J = aE, but instead the phenomenalogical 
equation of F. and H. London should be employed. Thus, 


E = iy3' 2 — 
p ° at 


( 6 ) 


The parameter (3“* is the London penetration depth, the distance through which the 
magnetic field falls to 1/e of its value at the surface of a superconducting half-space. 
It follows from equations (1) and (6) that 


d_ 

dt 


J(r) - J(0) - 0 



r 


H(r") • (a z x dr") 


= 0 


(7) 
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The quantity within the brackets in equation (7) is the well-known London fluxoid. Thus, 
equation (7) expresses the principle of fluxoid conservation. Starting from the London 
zero field initial conditions, it is seen that the fluxoid associated with this or any other 
contour within the superconducting cylinder is zero. This is the exact analog of the zero 
flux condition encountered in the case of perfect conductors. Conclusions 1 and 2 which 
were stated for perfect conductors remain the same for superconductors except that the 
word surface must be deleted. This is apparent since the infinite current density situa- 
tion associated with a time rate of change of flux no longer exists in this case. 

The magnetic field at a point r" in either cylinder can readily be found by summing 
over the contributions made by the current elements in both cylinders. Clearly then 


H(r”) 



( 8 ) 


But 


r” - r* 
|r*' - r» 1 2 


= + Vln r" - r'l 


where V operates on the field point variables. Therefore, 


H(r”) = 


©,(D 


2^ X ” ** I j d 2 r* 


When use is made of the vector identity 


V X SV = S(V X V) + VS X V 


it follows that 


H(r”)= / 


<D, ® 


J(r') v x In | F” - r* | j d 2 ? 


( 9 ) 


( 10 ) 


Since the integration is being performed over the source points, and V x operates only 
on field points, the curl operator can be taken out from under the integral sign. Hence 
in the cylinders, 
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(11) 


H(r M ) = vx - — 

L 2n 

It should be noted that equation (11) could also have been derived by first determining the 
vector potential for the system. If the vector potential A satisfies the requirements of 
the Coulomb gauge (V • A = 0), it can be shown that /u" 1 A is just the expression con- 
tained within the brackets in equation (11). Using the identity 

H(r") • (a z x dr") = ^(r’) ln l ?t ' 1^' dr " ( 12 ) 


and the fact that the fluxoid is zero yield 


L 


(D * 


J(r’) In r 


r* T - r* Id^r* 


J(r) = J(0) - £- 

2tt 


J(F') In I? Id 2 ? 




(13) 


Equation (13) is equivalent to that derived by Cooper and Marcus. The first two terms on 
the right-hand side are independent of r, therefore the current density can be expressed 
as 


J(?) = C +~jC 
2t r J(T\ 


2it © 


J(r’) In |r"- r* Id 2 ? 


(14) 


where C is an arbitrary constant determined by the total current. 

When the steps in this derivation are reconsidered, the following interesting fact is 
apparent: The nonuniformity of the current density in the film (J(r) ^ J(0), r" / 0) is due 
to the magnetic flux crossing the area bounded by the dotted contour shown in figure 1(b). 
For certain cross-sectional geometries, such as the strip transmission line which will 
be discussed later, this flux is small if the conductors carry equal but oppositely 
directed currents. Thus, in these cases, only small variations of the current density 
are expected. Furthermore, in the case of a strip transmission line, this concept 
allows certain properties of the kernel of the integral equation to be deduced so that a 
closed form solution to equation (14) can be demonstrated. 

Before proceeding to this solution, consider the case in which conductor 1 is super- 
conducting and conductor 2 is normally conducting. As was shown previously, the cur- 
rent density in the normal conductor is uniform (in the static case) so that if Ij ^ are 
the net currents in conductors 1 and 2, respectively, equation (13) becomes 
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( 15 ) 


J(r) - J(0) =~ / J(r') In |r - r* |d 2 r' / J(r') In |r’ |d 2 r» 

27r 4) 


+ 




r T 


d 2 r’ 




It is interesting to note that for this case the distribution of the current in conductor 1 
is not influenced by the presence of conductor 2 for I 2 = 0. This was not true for the 
case of two perfect conductors or two superconductors, as was shown previously. 


STRIP TRANSMISSION LINE OF RECTANGULAR CROSS SECTION 

Consider the strip transmission line, shown in figure 2, that consists of two parallel, 
infinitely long, superconducting thin films of rectangular cross section. Each is of 
width w and thickness d, and the centers of the two are separated by a distance i. 



s — W 

Figure 2 . - Strip transmission line. 


The parameters of the system are chosen to be in a range of practical interest charac- 
terized by 0. 0001 < d/w < 0. 01, and (/3d) <1. A current I flows into the top film, and 
an equal but oppositely directed current flows in the bottom film. For the case of the 
strip transmission line system just described, equation (14) becomes 


J(x,y) = C+f- 
477 


r w/2 rd/2 

cAw/2 J-d 


J(x’,y’) In | (x - x’) 2 + (y - y') 2 |dx' dy* 


-d/2 
’w/2 Md/2)-i 


rw/2 At 

f- / / J(x’,y') In |(x - x») 2 + (y - yO^dx’ dy’ 

•'-w/2 J-(d/2)-H 


(16a) 
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Using the change of variables x' — - (x’ + l) in the last integral on right-hand side in 
equation (16a) and noting that 


J(x,y) 


d . .d 

-<x<- 

2 2 


= - J(-x-£ , y) 


d d 
--<x<- 
2 2 


(16b) 


yield 


J(x,y) = C +y- 
4ll 



_[a -x-)V ly-j') 2 J.., dx , dy 

|(x + x’ + j 0 2 + (y-y') 2 [ 


(17) 


Since (/3d) < 1 it will be assumed that J does not vary appreciably in the x-direction in 
the film. (This assumption will be examined more closely in the next section. ) 

Therefore, J(x, y) « J(0,y) = J(y). It is convenient in equation (17) to perform the integra- 
tion over x' first. Let 


G(x’, y-y’) = fin |(x») 2 + (y - y’) 2 |dx* =x' In |(x’) 2 + (y - y’) 2 |-2x' + 2(y - y») tan" 1 ^^7 
Then, 


*d/2 


'-d/2 


] (x') 2 + (y - y’) 2 1 
ln . t ty — y i . i- dx t = G 

(x* + £) 2 + (y - y') 2 | 


- G 


d 

2 


d 

f+ 2 


2 


(18) 


Thus, 


f 

w/2 


'w/2 

J(y) = c + 4n I J (y’) d yM 


[d ln 1 4- + (y - y') 2 -2d + 4 (y - y') tan -1 1 - 
L I 4 y - y'J 


( 4 + 2) 111 1 ^ + + < y * y') 21 


is . + - \ 2 

- 2 (t +^)+ 2(y - r) ( £ " 2) 111 ( £ ' 2) + (y ‘ y,)2 + 2 ( £ ' 2)' 2(y ‘ y,) tan 


/I-4N 

tan'* ( — 

\ y-y) 


(19) 
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When the dimensionless variable u = — y is introduced, equation (19) becomes 


J(u) = C + 



J(u’)K(u, u’)du* 


(20a) 


where 


K(u, u’) = K( | u - u' | ) = In 



+ (u - u ’) 2 


+ 2 


u - u* 
d/w 


tan - -*- 




(20b) 


The kernel of the integral equation is an explicit function of the thickness to width 
ratio of the films, as well as the ratio of the separation distance to the width. It was 
mentioned earlier that the nonuniformity of the current density is due to the magnetic 
flux that crosses the area bounded by the dotted contour shown in figure 1(b). Certainly 
line currents (henceforth referred to as source points since the analysis is two dimen- 
sional) at all points in the cross section contribute to the total flux. In figure 3 let the 


H FP t 


( 

r 



■ Q i 


X P! 


O P 2 



Figure 3. - Flux contributions from antisymmetric source points. 


currents at Pj and Pg be equal in magnitude but oppositely directed. Let be the 
distance between points F and Pp and Q 2 be the distance between points P^ and P 2 . 
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If H-cvr, and H™, are the magnetic field components at F due to the currents at 
Fi l **2 

Pj and Pg, respectively, then if 


H 


FP, 


- H 


FP c 


H 


< 0.01 


FP. 


and 


0. 989 < cos 6 < 1. 000 

Thus, due to the dimensions of the transmission line that are being considered, contri- 
butions from antisymmetric source points (in the two films) at large distances from the 
observation point tend to cancel themselves out. It should be expected, therefore, that 
virtually the entire net flux will be contributed by source points that lie within some 
small distance from the observation point (x, y). This distance should be of the order of 
10 ft, where H is the separation distance between the two films. Thus the kernel should 
be sharply peaked about u = u* and, based on the preceding discussion, it should be ex- 
pected that the magnitude of the kernel will be a monotonically decreasing function of 
| u - u* | . Since SL >d and -(d/2) < x,x' < (d/2), the two-dimensional kernel in equa- 
tion (17) is negative definite. When the process is considered by which equation (20a) is 
derived from equation (17), it is apparent that the one-dimensional kernel in equa- 
tion (20b) is also negative definite. 

A plot of the absolute magnitude of the dimensionless kernel against (u - u 1 ), which 



- 95 Percent of total area 


Figure 4. - Absolute value of kernel of inhomogeneous integral equation (not to 
scale). 
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II 

1 


summarizes the conclusions of the previous discussion, is shown in figure 4. This ap- 
proximate sketch is supported by numerical analysis and shows that the kernel drops off 
to ~l/e of its maximum value when |u - u’ | is on the order of l /w, and that the bulk 
of the area under the curve (95 percent) is contained in the interval 0 < |u - u’|< io (£/ w) 
Note that approximately 50 percent of the area under the curve is contained in the range 
0 < |u - u* | < (f/w). 

If 



(21a) 


equation (20a) can be rewritten as 


J(u) = C - X 



(21b) 


Equation (21b) can be solved by the Liouville-Neumann method of successive approxima- 
tions. Let Jq(u) = C. Carrying through the integration on the right-hand side in equa- 
tion (21b) yields 


Jj(u) = C 


1 ~ x J f 1 | K d u - u il)| du i 


(22a) 


Repeating the process a second time but now inserting Jj(u) in the integrand of the 
right-hand member yields 


J 2 (u) = C 


1 - X 


f |K(|u- uJjJdUj 


+ X 


/i ii | K( I U 1 - u 2l)| | K <I U - u ll>| dugdu^ (22b) 


Thus, after n repetitions it can be shown that 
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1 


J n (u) = c 


n 111 

> + J ■/ | K <l u m-l - u m l>| - | K (I U - «il) 

m=l -1 -1 -1 


du m--- du l| 


(22c) 


Thus the current density is expressed as 

J(u) = lim J (u) (22d) 

n-oo 


provided the series in equation (22c) converges as n— °°. The range of convergence for 
this series will be discussed later. 

It is instructive to note that to zeroth order the current density Jq(u) is assumed 
to be uniform over the width of the film. The first-order correction to this assumption 
is included in J^u) and is computed by considering the sum of the interactions of all the 
source points in the films on a particular observation point. This sum 


X 



= X 



(23) 


where k(u) will be referred to as the coupling factor. The second-order correction, 
contained within ^(u), is determined by considering the effect of all the source points 
on a particular source point before obtaining the coupling factor. The higher-order 
terms in X are further expressions of the interactions of source points with source 
points. 

It is noted that the terms in the kernel are of essentially two types: 


In 



+ (u - u') 2 


and 


“ - Hi ^-1 a (d/w) 
d/w u - u’ 


(24a) 


(24b) 


The first term (expression (24a)) is easily integrated with respect to u', If 
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V = u - u' 


then 


/ in l a2 © t<u - u ' > i du,= / 

~ l 1 *^u-l 




dv = v In 




- 2v + 2af- N j tan -1 


a(d/w) 


u+1 


1 . 


= in l a - 2 ( d / w ) 2 t ( u + 1 )' 


I u+1 


u- 1 a 2 (d/w) 2 + (u - 1) 


lu-1 


4 + 2a [-Vtan” 1 


u + 1 


a(d/w) 


tan" 


u - 1 
a(d/w) 


= A, 


(24c) 


When integrating expression (24b) it is important to remember that the inverse tangents 
are restricted to their principal values; therefore, 
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r 


Thus 


/c(u) = A x 




a=2k-l 



+ 2A« 

a=2k+l 


a=l 


+ A 


2 


-A, 


a=2k-l 


a=2k+l 


(24e) 


where k = i/d. 

It is useful at this point to examine some of the properties of /c(u). Figure 5 shows 
the absolute magnitude of the kernel plotted against u* for various values of u. The 



area under each curve is |/c(u)| . Since the kernel is negative definite, k(u) is also 
negative definite. Due to the narrow effective width of the kernel (approximately 
10 (1/ w)), it is seen that k( u) is a weak function of u [(k(u) ~ k(0))] in the central region 
of the film. When u is within 10 (£/w) of the edge of the film, k(u) is a strong function 
of u, and in fact it can be shown that at the edges 


|/c(±l)[ |/c(0)| +e 

where € is a small positive number. These observations are summarized in figure 6. 
Physically this implies that since an observation point is only affected by those source 
points within a range «10 (£/w), the observation points which are not within 10 (£/w) of 
the edges are effectively in an infinitely wide film. 

When a composite of figures 5 and 6 is considered, it is readily seen that 



K(|u- uj) 


duj * 


«< 0 ) 


ku) 


(25a) 
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MO)|/(pd ) 2 " (sL 


Thus, in general 



u m l) '•* K (I U " u il) 


du 


m 


du i » 


K(0) 


|m-l| 


k(u) 


|*c(u)l 


(25b) 


k<0)l 



4 


-t 

i 

i 

i 

_L 



T^0.5MO)| 


1 u 



Figure 6. - Absolute value of coupling factor (not to scale). 


Equations (25a) and (25b) are very 
nearly true in the central portion of 
the film. The error in these approxi- 
mations, when u is within 10 (i/w) of 
the film edges, will be discussed in 
the next section. Combining equations 
(22) and (25) yields 


J(u) = C 



(-l) n [x|/c(0)|] 




(26) 


If x|/c(0)| < 1, the series in equa- 
tion (26) converges absolutely, and the 
current density distribution can be ex- 
pressed in closed form as 


J(u) = C 


M«(u)l 

1 + X | k(0) 


(27) 


Figure 7 indicates the range of values 
for the parameters d, w, and i over 
which equation (26) converges. For 
example, once a film has been de- 
posited /3, d, and w are fixed, and 
the range of convergence for equa- 
tion (26) with i/d is to be determined. 
Clearly, the value of i/d cannot be 
less than or equal to 1, for this would 
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Figure 8. - Current distribution (not to scale). 


From equation (27) and the relation 
netic field can be readily deduced: 


imply that the two films touch or overlap. Since 
Xj«(0)| < 1, the upper limit for SL/d can be 
found by setting x|/c(0)| = 1, dividing by (j3d) 2 , 
and finding that point on the appropriate d/w 
curve that has this value as its ordinate. The 
projection of this point on the SL/d axis yields 
the upper critical value for SL/d. Therefore 
1 < (SL/d) < (SL/ d)__. Following somewhat similar 
lines, one can fix any two of d, w, and i and 
determine the allowable values of the third 
parameter. If X |/c(0)| > 1, the series does not 
converge and the Liouville-Neumann method is 
not applicable. With the help of figure 6 an ap- 
proximate normalized sketch of J(u) against u 
can be drawn (fig. 8). 

H = /3~ 2 ( v x J), the x-component of the mag- 


H x (u) = 4 (V X J)x = 4 J(u) = 

x /3 2 x /3 2 c, y 


— - J(U) = 
/S 2 w du 


2CX 


/3 2 w(1 + a|k(0) | ) 



r rl 

_d 

/ |k(|u - u'|)|du' 

|) du 

\J_l \ 1 1 1 


(28) 


Removing the absolute value signs from the kernel, applying Leibniz's rule, and re- 
membering that 


■— f ( u “ u ’) = " ~ f ( u “ u ’) 
du du’ 


yield 


H x (u)=~^ ^7 r f dK (|u-u'|) 

P w(l + x|/c(0)|) J- 1 

H x ( u ) = 

P wi 


K(|u- 1|). 

K( 

u + l|) 

1 + X 

K(0) 



(29) 

(30) 
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where the plus and minus signs refer 
to the top and bottom films, respec- 
tively. An approximate sketch of 
equation (30) is given in figure 9. 

In principle, Ampere’s law can be 
used to evaluate H y (x, u); however, for 
the purposes of the discussion in the 
next section, only the form of H y is 
necessary. In the central portion of 
the films it can readily be seen that 
H is maximum at x = -(d/2), 

-p - (d/2)J and approximately zero at 
x = (d/2), -[£ + (d/2)]. 

ERROR ANALYSIS AND CONCLUSIONS 

There are two sources of error in the analysis presented in the last section. The 
first error was introduced when it was assumed that the current density does not vary in 
the x-direction for the case (|3d) < 1 (one-dimensional approximation). The second error 
is associated with the approximations in equations (25a) and (25b) (the edge approxima- 
tions). 




It*. 


r L^ [ 
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* Mi ITT 1 


!4 


II 1 " 


G I 


7 r^L~ 


Figure 9. - Vertical component of magnetic field (not to scale). 


One-Dimensional Approximation 

In order to evaluate this error it is convenient to use a self-consistency argument. 

In other words, the one -dimensional solution (eq. (27)) is resubstituted into the two- 

dimensional integral equation (eq. (17)) in order to find the x-direction variation of J. 

In practice, it is more convenient to use the differential counterpart of equation (17), 

2 2 

namely V J = /3 J. The x-direction variation is then used to find an improved 
y-direction variation. For the purposes of this report it is not important if J is not 
uniform in the x-direction, as long as the inclusion of this variation does not appreciably 
affect the y-direction variation. It is shown below that this is, in fact, the case. 

From the behavior of the magnetic field in a superconducting film it can be argued 
that the largest x-direction variation of J occurs in the central portion of the film. 

This then is the worst region as far as an error in the one-dimensional approximation is 
concerned. Since both the current and the magnetic field satisfy linear differential 
equations of the same form, it follows that in the central portion of the films H y (x) must 
be a solution of the equation 
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(31a) 


d 2 H(x) 

—V ■ *y« 

dx 


subject to the boundary conditions mentioned in the last section: 


H y (x = + d/2) « 0 
H y (x = -d/2)-J w 


(31b) 


It is easy then to show that in the central region of the film 


w 


H y (X) sinh(/3d) 


sinh /3|x - (d/2)J 


(31c) 


Consequently, 


j m ) = ^ . j(o> g ? sh % a/ 2 ) ] 

z dx cosh /3(d/2) 


(32) 


It is seen that for (/3d) = 1, ~ 1. 5 J^j. Thus there can be an appreciable 

x-direction variation. 

Substituting equation (32) into equation (17) leads to an improved two-dimensional 
integral equation: 


/>w/2 /*d/2 

J(x, y) = C + — / J(y')dy' / 1 ln J{ x - x'> 2 1 (y - yjfl 

4 "4/2 id/2 cosh <W/2) (x + x' + J ) 2 + (y - y', : 


dx’ 


(33) 

To simplify comparison with the results obtained previously where J was assumed con- 
stant in the x-direction, it is advantageous to set x equal to zero in the integration over 
x*. Furthermore, it is clear that the largest contributions to the kernel come from 
source points closest to the field points. Thus in this ,, worst ,, region it is helpful to set 
y = y’ in an attempt to establish the upper and lower bounds for 
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cosh /3[x* - (d/2)] 
cosh(/3d/2) 


x 


,2 


(£ + x') 2 


dx' 


(34) 


Since every function can be expressed as the sum of its even and odd parts, expres- 
sion (34) can be rewritten as 



cosh(/3x’) In 



(-1) tanh 



sinh(/3x’) 


In 



dx’ 


(35) 


Examination of the integrand indicates that the second integral in expression (35) is 
always positive. If 


f( x t\ _ cosh /3 [x 1 - (d/2) 
cosh(/3d/2) 


and 


x« 2 

g(x') = In 

(fc + x') 2 

f d/2 

then the first integral in expression (35) is just y f E (x , )g E (x’)dx’, where f E (x') and 

g E (x’) are the even parts of f(x') and g(x’), respectively. Similarly, f 0 (x’) and g Q (x f ) 
are the odd parts of f(x») and g(x'). If is the average value of f E (x’) in the 

interval - - < x' < then 
2 “2 

f E (x ’> = % + m E (x,) 

Thus, expression (35) can be rewritten as 



f(x , )g(x’)dx' 



m E (x , )g E (x')dx’ + positive term 


( 36 ) 
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Examining each of the integrands reveals that 


/d/2 

/d/2 


f(x , )g(x')dx t < 0 


/d/2 

/d/2 


g E (x , )dx' < 0 


/d/2 

/d/2 


m E (x , )g E (x , )dx' > 0 


(37a) 


(37b) 


(37c) 


Thus, 



f(x*)g(x»)dx» > f E 



(38) 


Since 


/•d/2 

/d/2 g o <X ' )dX ' = ° 


and g(x f ) is always < 0, expression (38) yields 
•d/2 


f(x»)|g(x’)|dx» < f E jf d/2 |g(x')|dx' = S -- ^" /2) y^ 2 |g(x’)|dx’ 


(39) 


In a similar manner the lower bound can be found for 


JJ f(x , )|g(x , )|dx f 


-d/2 

The result is equal to expression (39) with |sinh(/3d/2)/(|3d/2)J replaced by unity. Since 
/d/2 

/ g(x/x , ,y/y , )dx’ can be identified with the kernel of the integral equation (eq. (19)), 

/ d/2 

it follows that use of the self-consistency procedure yields an improved u-variation 
integral equation 
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where 


(40a) 


J(u) = C - aX 




1 < a < sinh /3d/2 
0d/2 


(40b) 


Thus, at worst, a = 1. 04 (for the case /3d = 1) in the central region of the film. If 
X is replaced by ax in the closed form solution, the magnitude of J(u) in the central 
region is less than 2 percent smaller than that predicted by equation (27). As was 
pointed out previously, this error is less near the film edges. 


Edge Approximation 

It would appear from figure 6 that while the approximation /c(u) » /c(0) is a reason- 
able one to make in the central region of the film, its validity, as well as its usefulness 
at the edges, is questionable. Though the closed-form solution might therefore not be 
strictly true near u = ±1, equation (27) could be used as a reasonable trial function in an 
attempt to obtain a numerical solution for the one-dimensional integral equation. 

Numerical solutions were obtained in this manner with the aid of an IBM 7094II 
computer. The results revealed that the greatest error in the closed form solution oc- 
curred for k k c , where J(±l) was approximately 10 percent smaller than the value 
given by equation (27) for films with the same parameter values. It should be noted 
however that for cases where k < 0. 65 k c , this discrepancy is reduced to less than 
5 percent. 

If it is assumed that the superconducting films switch to the normally conducting 
state when J at any point in the films exceeds a critical value J c , it is apparent that 
the switching is initiated at the film edges. From the error analysis it is clear that if 
the J(u) given in equation (27) is used to calculate a critical current I (in terms of J c ), 
this value for the critical current will be smaller than the true critical current by less 
than 10 percent. It is noted that in calculating I the x-direction variation of J in the 
central portion of the film should be taken into account. This is easily done by combin- 
ing equations (27) and (32) as is done in the sample critical current calculation shown in 
the appendix. The critical current density can be determined by making use of a micro- 
scopic theory (ref. 7) and/or by experimentally determining I c for a particular choice 
of parameters which lie in the range treated in this report. Once J c is known, I c can 
be determined by using the above calculations for any other set of parameters which lie 
in the range under consideration in this report. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 6, 1965. 
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APPENDIX - SAMPLE CRITICAL CURRENT CALCULATION 

The critical current I of a superconducting strip transmission line, for the case 

V 

d/w « 1, can easily be calculated. For this case, the current is uniform over a very 
large percentage of the film width so that from equation (32) it can be shown that 

I ~ J(0) w(/3" X ) 2 sinh(/3d/2) (Al) 

The dimensionless variable u will not be used in this section. From equation (27), and 
making use of the fact that Jk(±w/2)| » ~ |/c(0)| , it can be seen that for a|/c(0) | < 1 


J(0) ” 1 + M k<0) I < a2 > 

where the variation of J in the x-direction at the film edges has been neglected in line 
with the discussion in the last section. According to the critical current density hypo- 
thesis, I = I„ where J(±w/2) = J„. Using this hypothesis and combining equations (Al) 
c c 

and (A2) yield 



J c wO -1 )2 sinhOd/2) 

1 + 2 ^ I K (0) I 


(A3) 


The numerator of the right-hand side of equation (A3) is the critical current I calcu- 
lated on the assumption that the current density is uniform in the y-direction. It can be 
seen from equation (24) that 


From figure 8 it is evident that equation (A4) is satisfied for the case d/w < 0. 001. It 
is interesting to note that in this range A k ( 0 ) is independent of (d/w). Therefore, 


cu 



(A 5) 
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For situations in which (d/w) is not much less than unity, equations (Al) and (A4) are not 
valid. In these cases the general expressions describing the variation of current density 
in the films can be used to calculate the critical current. 
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